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Abstract 

An expression for the joint probability distribution of the principal curvatures at an arbitrary point in the 
ensemble of isosurfaces defined on isotropic Gaussian random fields on W is derived. The result is ob- 
tained by deriving symmetry properties of the ensemble of second derivative matrices of isotropic Gaussian 
random fields akin to those of the Gaussian orthogonal ensemble. 
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I. INTRODUCTION 



We closely follow the notation in [ ] and [2]. Let T be a tensor with rank a and dimensions 
(di,..., da). The bijective linear map vec associates T to the vector vec T in R N , N = YYLi 
with entry (vec T) k , k e { 1 , . . . , flti d i)> § iven b Y ( vec T )t = T *i wim ij e {1, . . . , dj\ uniquely 
defined by k = 1 + Z"=i( z ; _ 1) Y\iZ\ dk- Let d t = n, i = 1, . . . , a. The linear operator diag maps 
T to the vector diagT in R" with entry (diagT)^, k e {1, . . . ,n}, given by T kk j,. If restricted to 
the set D of diagonal matrices, D 3 D — » diag D is bijective, and therefore the inverse mapping 
diag 1 is well-defined. Let a = 2 and d\ = d 2 = n. The linear operator uni maps T to the vector 
uniT e R"(«+i)/2 w ^jj enrr y (vecT)*, k e {!,..., n(n + l)/2}, given by (vecT)* = Ty uniquely 
defined by = (_/ - l)n - j(j - l)/2 + i. This operator maps a matrix T to a vector containing the 
entries of T read along its columns, but ignoring the elements above the main diagonal. 

Following the notation in [ 1 ], the (n, n) identity matrix is denoted by I„, C„ denotes the (n 2 , n 2 ) 
commutation matrix, i.e., the unique (n 2 ,n 2 ) matrix such that vecT T = C„vecT for all (n, n) 
matrices T. The matrix D„ denotes the (n 2 , n(n + l)/2) duplication matrix, i.e., the unique (n 2 , n(n + 
l)/2) matrix such that vec T = D„ uni T for all (n, n) symmetric matrices T. From this definition, 
we have uni T = vec T, where A + indicates the Moore-Penrose inverse of the matrix A. The 
duplication matrix and the commutation matrix are related through the identity D„D+ = (1/2)(I„2 + 
C„). Finally, l, v „ denotes an (n, m) matrix of ones. 

A random scalar field in the set X is a stochastic process, i.e., an indexed set Tx = 6 X} 

of random variables T x defined over the same probability space (f2, cr n , P). A random tensor fields 
is defined analogously, with T x a vector- valued function. 

Let T x and T x be two random scalar fields as above, and assume that is zero mean, which, 
for the purposes of this work, implies no loss of generality. If the expectation E{T^T 2 } taken 
over the joint distribution of T\ and T% is defined for all x and y in X, the function R r , r2 (x, y) = 
E{T^T 2 } defines the cross-covariance function of the random fields. If ^ = T x = Tx, the 
notation is simplified to Rf — Rr,r, an d the function Rf is referred to as the autocovariance 
function of Tx- The conditional autocovariance function of T x given T%, Rf^r 2 , is defined by 
taking the expectation of T x l over the conditional distribution of T^ x and Tj given T%. In the case 
of a random tensor field the definitions are analogous, with the product T^T 2 replaced by a tensor 
product vec T x l <8> (vec Ty ) T and the expectation taken over each entry of the tensor. Henceforth, 
the term "random field" will be used in reference to both scalar and tensorial random fields. 
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Let X be a vector space. Zero-mean random fields T x and onl are wide-sense stationary 
if their cross-covariance function R r i r 2(x,y) satisfies ^(x,y) = Rr x (s) with s = x - y for all 
x and y in X. Henceforth we will assume that X = W, and therefore the subscript X in Tx can 
then be safely omitted by defining T - Tw 1 - A stationary random field on R" is isotropic if its 
autocovariance function Rf{s) satisfies Rf{s) = cr 2 p(||s||), where p is a correlation function [2] and 
| HI is the standard Euclidean norm in R". 

H. DERIVATIVES OF ISOTROPIC RANDOM FIELDS 

Great simplification is achieved in the derivations that follow if p(||s||) can be rewritten as 
p(||s||) = r(||s|| 2 )/cr 2 <=> p( VlMI) = Klls||)/cr 2 . For ||s|| > the smoothness of r(||s||) is contingent 
upon that of p(||s||). However, for ||s|| = the non-differentiability of VlMI could be a problem. 
This is not the case, as shown in appendix. The symbols p® and r[, denotes the z'-th derivative of 
p(||s||) and r(||s||) with respect to s = ||s|| at s = 0. 

Let / : R — » R be a scalar function. The symbol is used to describe the matrix of 

dimensions (n a , n b ) of partial derivatives of /, i.e, 

/ d a+ V \ d" +b f 

\ dx" T dx h ) A B dx ai . . . dx a dx h . . . dx h ' 

where a t e {1,2, ... ,n] and bj e {1, 2, . . .,n] are uniquely defined by A = 1 + Yli=\( a i ~ l)^'" 1 
and 5=1 + 2" =1 (/?y - l)n ;_1 . Second differentiability of the autocovariance function R r (x,y) 
of a random field at the pair (x, x) implies mean square differentiability of the random field itself, 
as demonstrated in theorem 2.4 of [2]. If a mean-square differentiable random field is stationary, 
its derivatives will also be stationary. For T l and T 2 jointly stationary with cross-covariance 
R r i >T 2(s), we define R° rl f2 = R rl>r 2(0). 

The theorem that follows is central to this work: 

Theorem 1. Let T be an isotropic random, field on R" with autocovariance function Rf(s) = 
cr 2 p(||s||) where p(s) is four-differentiable. Let dT and d 2 T be the tensor fields defined as 

4 or „ 2 A o 2 r 

dr =— and 5f=— . 
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Then 



4r = -vVoX, (la) 

R% r = ^ 2 Po 4) (I„2 + C » + vec J " vecT l »>- (lc) 

Proof. Using lemma 4 in appendix we write p(||s||) = r(|]s|| 2 )/cr 2 . From the four-differentiability 
of r and theorem 2.4 in [2] we have 

d 2 r(Hs|| 2 ) 

flsrO) = — d , (2a) 
<9 3 r(||s|| 2 ) 

Rar&ris) = ~ dsTgs2 . (2b) 
d 4 r(||s|| 2 ) 

= • (2c) 

The chain rule can be used to expand (2a)-(2c), and substituting the identities || = ^ = I„, 



= I„ ® I„ and = C„ in the result produces 



Rsris) = -4r (2) (||s|| 2 )s ® s T - 2r (1) (||s|| 2 )I„, 
#ar,#r(s) = -4{2r (3) (||s|| 2 )s ® s T ® s T + 
r (2) (||s|| 2 )(s T ®I„+ 
I„ ® s T + s ®vec T I„)}, 
R a i T (s) = 4{8r (4) (||s|| 2 )s ® s ® s T ® s T + 
8r (3) (||s|| 2 )(vecL,®s T ®s T + 
2s <g> (I„ ® s T + s T ® I„) + s ® s ® vec l£)+ 
r (2) (||s|| 2 )(I„ g) I„ + C + vec I„ ® vec T !„)}. 



Making s = completes the proof. 



III. CURVATURES OF GAUSSIAN RANDOM FIELDS 

Let /(x) be a second-differentiable scalar function on R". For each x € M." for which df/dx £ 
we define the set F x as F x = {y e R" such that /(y) = /(x) and d//dy * 0}. If F x * 0, F x 
is a (n - l)-hypersurface in R" [3]. Let df = df/dx T and <9 2 / = d 2 f/dx T dx. The principal 
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curvatures of the hypersurface F x at x are given by the set of eigenvalues k obtained by solving 
the eigenproblem 

/ df®{dff \ d 2 f 

L ^r-r:^ \-rrr-r: V = KV , (4) 



\W\\ 2 ) \\df\\ 

with v e R", ||v|| = 1 and \ T df = [4, pg. 138]. Let {n,-, i = 1, . . . ,n - 1} be an orthornormal basis 
for the null-space of df, i.e., njdf = and njiij = Sjj, and let N be the matrix N = [ni . . . n„_i]. 
The eigenproblem in (4) can be rewritten as 

N T (d 2 /)N 



U = KU, (5) 



x 



with u e W~ l , ||u|| = 1. Equation (5) is still valid if the function / is the realization T(to), 
co e Q, of a mean-square second-differentiable scalar random field T on W. Therefore the random 
tensor field 'K of curvatures of isotropic Gaussian random fields is implicitly defined at x such that 
&F x (aS) + by the solutions of the equation 

N T (d 2 T)N 

where N is a random tensor field satisfying NjdTx = 0, and NjN x = I„_i, and 1A is the tensor 
field such that U x (co) are the eigenvectors associated to the eigenvalues ( K x {co). 

Henceforth we assume that T is an isotropic, second-differentiable Gaussian random field, 
which is defined simply as an isotropic random field for which the joint distribution of any finite 
set of random variables {Ta), A 6 K", is Gaussian. This assumption implies that the zero-mean 
random tensors dT and d 2 T are also Gaussian, and therefore &F X and d 2 T x are fully characterized 
by their covariance matrices, given by (la) and (lc) in theorem 1. However, because d 2 ^ is 
symmetric, its probability density must be handled with care, since R® 2r is not invertible. 

The following lemma is a trivial corollary of the theorems in [1, sec. 7]. 

Lemma 1. Let A„ be a (n, n) invertible matrix. Then 

(A* ® A n )R% r (A; 1 ® A,; 1 ) = R% r . (7) 

Let d 2 T be as in theorem 1, and let K be a (n, m), n > m, orthonormal tensor field independent 
of d 2( F, i.e., a random tensor field such that for all x 6 R' ! any realization K{co) of K satisfies 
<RI(uj)<R x (u)) = I m and <R X is independent of d 2 T x . Define d 2 T' = < R J d 2 r , R = {n x d 2 T x K x , xel"). 
We prove the following lemma: 



5 



Lemma 2. d 2 T' is a Gaussian random tensor field with autocovariance function Rgiyr,(s) such 
that 



R% r , = o- 2 p { * ] (I m i + C m + vec l m vec 1 I m ). (8) 



Proof. Let P d i r ^ Kx be the conditional probability of the random tensor d T x given % x . Since d T x 
is independent of ^ x , d 2 T^ given K x is a linear function of d 2( F x , and therefore it is zero-mean 
Gaussian. Using the identity vec(A5C) = (C T <8> A) vec B and properties of commutator matrices 
[1] we can write vec d 2 T x = (R X T 9 K T ) vec d 2 T x , and therefore 



R% niK% = (K T ® K T )R% r (K 9 (9) 
= cr 2 po 4) (I m 2 + C m + vec l m vec T I m ), (10) 

using lemma 1. Since d 2 T x given K x is zero-mean Gaussian and R® 2r ,^ does not depend on 
^ x for fixed m and n, we have P d i r ^^ = P Therefore Pd 2 r; is zero-mean Gaussian with 
autocovariance function satisfying (8). □ 

Lemma 2 justifies the notation d 2 T m = < R T d 2 T'R, for % (n, m). Let 6 2 T n be the (n(n + l)/2, 1) 
vector defined as 6 2< F n - uni d 2 T n = D+ vec <9 2< F". Its covariance matrix H n is invertible and given 
by 

l„ = v:r% t ,x>^. (id 

Therefore the probability density pgi<j^ of d 2 !F x is standard: 

1 / h^hx 

/Wh) = —== exp -2— . (12) 



We now define the random fields (^",X") = {(^,£ x ) e 50(n) x R'\ x e R"| K' T ^ X = I» and 
^ x T diag _1 X x ^ x = d 2 T x ), where 50(n) is the special orthogonal group of (n,ri) matrices. Let 
eig" 1 be the mapping given by 

eig" 1 : SO(n) x R" -> S(n) 

(13) 

(R,A)^S = eig,; 1 (R,A) 

where S(«) is the set of (n,n) symmetric matrices. This mapping is differentiable and onto, and 
therefore the joint probability density p-R x x x °f and X x is given by 

p K>£ n(R, A) = /7^(uni(R T diag" 1 AR))\J(R, A)\, (14) 
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\J(R, A)\ is the absolute value of the Jacobian determinant of eig" 1 . 

Theorem 3.3.1 in [5] provides a "closed-form" expression of the probability density of the 
eigenvalues of random matrices in the Gaussian orthogonal ensemble (GOE„). This is the ensem- 
ble of (n, n) real symmetric matrices M. with probability density invariant with respect to similarity 
transformations h\ — » R T A\R for any given (n, n) orthonormal R and such that the probability dis- 
tribution of distinct entries are independent from each other. Even though each realization of d 2 T n 
is real, symmetric, and, by applying lemma 2, invariant to the same similarity transformations, 
d 2 T n is different from GOE„, because the distinct entries of &F£ are not independent. However, 
the assumption of independency used in [5] is important only to derive an expression for the joint 
probability density of the entries of random matrices in GOE„, and we already have that for the 
matrices in 5 2 T n . Once such an expression is available the result in [5] derives from the observa- 
tion that, in an expression analogous to (14), the term R appeared only on |/(R, A)\, and therefore 
the probability density of the eigenvalues of matrices in GOE„ is obtained through the integration 
of |/(R, A) | over SO{n). The next lemma shows that this is also the case for the probability density 
P£n of the eigenvalues of d 2 T£ : 

Lemma 3. Let A e W, R e SO(n), and £„ as in (11). Then 




(15) 



where 



t n = a 2 pf{2\ n + l^iy. 



Proof. The following identity can be easily verified: 




(16) 



Let A = diag 1 A and u R = uni(R T AR). Therefore, 



u/S^Ur = [uni(R T AR)] T E; 1 [uni(R T AR)] 

= [D; vec(R T AR)] T i:,; 1 \D + n vec(R T AR)] 
= (vec A) T (R ® R)D+' T E,; 1 D,;(R T O R T ) vec A, 
and, using [(R ® R)D„] + = D+(R T ® R T ), 

u/E;> s = (vec A) T [(R ® R)D„2:„D^(R T ® R T )] + vec A, 
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which, using D„£„Dj = R° 2r „, yields 

u/S^Uk = (vec A) T [(R ® R)#° 2r „(R T ® R T )] + vec A 
= (vecA) T (^ 2r „) + vecA, 

= (17) 
using lemma 1 and (16). □ 

The integration of |/(R, X)\ over SO(n), carried out in [5], gives 

/7r (H+D/4v« n"", 1 n" . , \a -a \ 



r 

J so 



'50(n) 

and it can be shown that the determinant of E„ is given by 

|E„| = 2"- 1 (2 + n)(cr 2 p! ) 4) r ( " +1)/2 . (19) 
Together with lemma 3, these results demonstrate the following theorem: 
Theorem 2. The probability distribution of the eigenvalues ofd 2 T£ is 

2(2-7n-« 2 )/4 

P£ ' M) = V2T^(o-^r ( " +1)/4 n;L 1 r ( i + //2) X 

n— 1 « 



j^PJ|A y -Ai|«p(-^|^). (20) 

i=l j=i+l 



Since N x in (6) is a function of &F£, theorem l(lb) implies that d 2( F£ is independent of N x 
for all x. Therefore dT^~ x = N x dT£N x according to lemma 2. Theorem 2 can then be applied to 
obtain an expression for the probability density of the eigenvalues of the numerator of (6), pj?-\. 
Using theorem l(la), we can show that the denominator of (6), ||<97^'||, is distributed according to 
cr(-p ) l l 2 X{ri), where Xin) follows a ^-distribution with n degrees of freedom, and therefore its 
probability density p\\ar£\\ is given by 

2m"- 1 exp[ a 2 /(2o- 2 pi 2) )] 
(-2^p (2) ) n/2 Y(n/2) 

We can now prove our main result: 
Theorem 3. Let 'K be as in (6). Then 

2 ( " 2 - 7 " +8)/4 r[n(n + l)/4] 

p%M) = . — : x 

VTT^r(n/2)n;Ci 1 r(i+//2) 



a 



n(«-l)/4 FT"- 2 FT" -1 



MKC/^-irS-r/^ + i}^ 

where a = -p| ) 2) /(2p[ ) 4) ). 



(22) 



Proof. Since &F£ and d 2 T£ 1 are independent, so will be ||<97\"|| and X* 1 ■ Using (6), we have. 

PkM)= \ u n ~ l p Ll -i{Ku)p\ W£ \\{u)du. (23) 
Jo 

Substituting (20) and (21) in (23), we obtain (22). □ 
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APPENDIX: DIFFERENTIABILITY OF THE AUTOCORRELATION FUNCTION 



Correlation functions are characterized by the Wiener- Khintchine theorem, a simplified version 
off which, shown below, is quoted verbatim from [2]: 

Theorem 4. A real function p(||s||) on R" is a correlation function if and only if it can be repre- 
sented in the form 

PWD - 2«-*IW2) £ |^jggW), (A.l) 
where the function <D(fc) on R has the properties of a distribution function and J v is a Bessel 
functions of the first kind and order v. 

Lemma 4. Let the i-th moment of the distribution O(fc) in theorem 4 be defined. Then, the i-th 
derivative ofr(\\s\\), r (i) (||s||), exists and is given by 

f«>(INI) = 2™T(d/2) J~ k iJ ~!0^- dmi (A.2) 

Proof. Define the operator D ; acting on a function f{u) as 

Di[f(u)] = (-4-) [/GO] (A.3) 



K udu i 

where the term in the right-hand side is recursively defined as 



I*) um = itm (A . 4) 

udu L u du 



--H [/(«)] = --r 
udu }■ \udu 



-T-l 

u dul j _ l 



(A.5) 
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It can be shown by induction that 



r (0 (||s|| 



1 d , 

--H lp(«)] 

u du 



(A.6) 



The operator D, and the integral in theorem 4 can be interchanged, since the functions and the 
measure d<&(k) involved satisfy the conditions of Lebesgue's dominated convergence theorem. 
The identity 



Di 



Jviu) 



(-D 



Jy +i (u) 



found in [ ], completes the proof. 



(A.7) 
□ 
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